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We calculate 0(a%m ) and 0(a%m log as) corrections to the quarkonium IS energy level an¬ 
alytically, which have been overlooked in recent studies. These are the contributions from one 
Coulomb-gluon exchange with the 1-loop vacuum polarization insertion. A part of the corrections 
was computed numerically some time ago; we correct its error. 
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Heavy quarkonium systems, such as bottomonium and 
charmonium, provide an important testing ground for 
studies on physics and dynamics of QCD boundstates us¬ 
ing non-relativistic boundstate theory based on pertur¬ 
bative QCD. Recently accuracy of the theoretical pre¬ 
dictions has improved dramatically. This owes to the 
progress in computations of higher-order perturbative 
corrections and the discovery J^] of cancellation of leading 
renormalons in these systems. So far, impacts in related 
fields have been as follows. (1) The mass of the bottom 
quark in the MS scheme has been determined accurately 
from the T spectrum |Q . It has various important appli¬ 
cations in B physics. Also it posed tight constraints on 
some models beyond the Standard Model such as super- 
symmetric Grand Unified Theories (2) The bottomo¬ 
nium spectrum is reproduced reasonably well by pertur¬ 
bative QCD up to the n = 3 levels. A new physical 
picture on the composition of the bottomonium masses 
has been proposed on the basis of the calculation [Q. 
(3) After incorporating the renormalon cancellation, the 
perturbative static QCD potential is shown to agree well 
with potentials of phenomenological models and with the 
QCD potential computed by lattice simulations ||]. 

Despite of these developments, in the context of non- 
relativistic boundstate theory, there have been no sys¬ 
tematic methods for identifying higher-order corrections 
to physical quantities of heavy quarkonia (energy levels, 
decay rates, etc.) in the non-relativistic expansion in 1/c 
(c is the speed of light). Rather, experts have identified 
separate contributions from inspections. A first attempt 
to compute 0(l/c 2 ) corrections to the quarkonium en¬ 
ergy levels was given in j(| , although its main results con¬ 
tain errors. |^| In terms of the (effective) non-relativistic 
Hamiltonian for quarkonium systems, different terms up 
to 0(l/c 2 ) have been identified as follows. The part 
which is identical with that for QED boundstates has 
been known for a long time (the Breit Hamiltonian). 
The 1-loop correction to the static QCD potential was 
computed in |?J; the 2-loop correction was calculated in 
H. The 1 /r 2 -type potential unique to non-abelian gauge 
theory was derived in |5| |9j. The analytic expression of 
the quarkonium energy levels have been derived from the 
above Hamiltonian in jIo| . It has been considered that at 
this stage all the corrections up to 0(l/c 2 ) = O(agm) to 


the energy levels have been identified and calculated, (as 
is the coupling constant of QCD.) However, the O(agm) 
correction from the non-instantaneous part of the gluon 
vacuum polarization computed in [Q has not been taken 
into account in all the studies which appeared after that 
computation. It seems that this part of the results of j(]] 
has never been taken seriously and has been forgotten 
largely. For instance, 0 states incorrectly that con¬ 
tributions from the non-instantaneous part of the gluon 
vacuum polarization are at most O(a^m). 

Recently many efforts have been devoted for identify¬ 
ing efficiently and systematically higher-order corrections 
in non-relativistic boundstate theory. For this purpose, 
people have constructed effective theories such as non- 
relativistic QCD (NRQCD) [jTTJ and potential-NRQCD 
(pNRQCD) Q, or used the asymptotic expansion of 
Feynman diagrams in 1/c (threshold expansion tech¬ 
nique) Ju|. As for the quarkonium energy levels, these 
formalisms provided more solid bases to the consensus 
that all the 0(1/c 2 ) corrections are taken into account. 
People have begun calculations of some of the corrections 
at 0(l/c 3 ) and beyond 0 Hi- There seems to be 
a tendency in the held that (in particular) the pNRQCD 
formalism is considered to be already well established for 
organizing and identifying higher-order corrections sys¬ 
tematically. Nevertheless, the 0(l/c 2 ) correction men¬ 
tioned above (and given in a more complete form in this 
paper) seems to be overlooked in this formalism, at least 
in its present usage. To our knowledge, up to date there 
exists no method which enables identification of all the 
higher-order corrections systematically. 

In we found a correction to the 15 energy level at 
0(agm log as) (within the approximation of a bubble- 
chain resummation). The correction originates from the 
ultra-soft region ko, \k\ ~ a 2 s m of the momentum of 
the gluon exchanged between quark ( Q ) and antiquark 
( Q ). From further inspections, we find corrections to 
the energy level at 0(agm log as) and 0(a A s m) origi¬ 
nating from the ultra-soft and the soft momentum region 
/c 0 ) |fc| ~ asm. Both of these corrections stem from the 
non-instantaneous nature of the gluon vacuum polariza¬ 
tion. In this paper we compute these corrections to the 
energy level of the quarkonium 15 state analytically in 
Coulomb gauge (see Fig. [|). We assume that Q and Q 
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have equal masses. The gluonic contribution corresponds 
to the correction A se Mig computed numerically in 
Following ]n|, our computation of the perturbative 
corrections is based on the Bethe-Salpeter (BS) formal¬ 
ism [|n| . We work in Coulomb gauge, since severe gauge 
cancellations occur in other gauges [|o). (Most com¬ 
putations in non-relativistic boundstate problems have 
been carried out in Coulomb gauge conventionally, ex¬ 
cept where gauge-independent contributions have been 
computed.) According to the BS formalism, the poten¬ 
tial energy E pot and the self-energy contributions E$e in 
the total energy of the boundstate can be expressed, in 
the rest frame of the boundstate, as 

E pot = (x- k -X) , ( 1 ) 

Ese = [(Sf,q)~ 1 (S Fi q)~ 1 ] ■ x) ■ ( 2 ) 

Here, XiX and E denote the BS wave functions and the 
BS-kernel, respectively; M is the boundstate mass (to¬ 
tal energy); Sf represents the full propagator of Q or Q. 
We set the momentum of the center of gravity as (M, 0). 
The dot (•) represents contraction of spinor indices and 
an integral over the relative momentum between <5 and 
Q. Diagranrmatically, E po t (Ese) represents the contri¬ 
butions from the diagrams where the Q and Q lines are 
connected (disconnected), and E pat + E$e = M. 

We compute the contribution to the potential energy 
from the 1-loop vacuum polarization of the Coulomb 
gluon n(fc): 

<? = (n(‘)) is 

i I" d A p f d A k 

~ 2 M[ 0 J J ( 27r ) 4 J ( 27r ) 4 

x (p) Kc(k) n(fc) x <°> (p + k). (3) 


Q 
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FIG. 1: The correction to the quarkonium potential energy 
induced by the 1-loop vacuum polarization insertion to the 
Coulomb gluon propagator. 


Here, m denotes the pole mass of Q and Q. <pc, i s (p) 
and m[°J denote the IS Coulomb wave function and its 
energy level, respectively: 


<t>c,is(p) = 7 


\f 7 ht {CfchsiE) 2 


p 2 + ( CFasm/2y 


2 ’ 


M[ 0 J = 2m- 


( Cfcxs ) 


A represents the spin of the boundstate, 

l + 7 o i_ 7 o 
A j = —-— <Jj —-—, 


(7) 

( 8 ) 

(9) 


where aj = 75 / 1/2 for J = 0 and a.j = //i/2 for J = 1 
and polarization vector e. 

First we compute n(fc) using dimensional regulariza¬ 
tion and in the MS renormalization scheme. (We assume 
that spacetime contains one time dimension and 3 — 2e 
space dimensions.) The Feynman diagrams contributing 
to n(fc) are shown in Fig. g. Since this is a 1-loop calcula¬ 
tion, it can be carried out more or less straightforwardly. 
Some useful techniques for calculations in Coulo mb g auge 
with dimensional regularization can be found in |2l| . We 
write 


This is shown diagrammatically in Fig. g. We assume 
that the quarks in the vacuum polarization are mass- 

(1 C\ 

less and compute SE^ ot ' in a series expansion in as 
up to 0(a A s m ) and O(agmlogas)- In Eq. (g), K C = 
—iCf ( 47 ras/|fc| 2 )( 7 °(g) 7 0 ) is th e kernel of one Coulomb- 
gluon exchange at tree level. Throughout this paper 
as = cts{p) denotes the strong coupling constant de¬ 
fined in the MS scheme with m active flavors, where p 
is the renormalization scale; Cf is the Casinrir opera¬ 
tor of the fundamental representation of the color SU (3) 
group. The BS wave functions in the leading-order of 
1 /c-expansion are given by 


X$(P) = TO + TO)] \/2M[s<PcMp) A, (4) 
Xisip) = [E>(p)+T>(-p)] \J2 m[° s ] (5) 


V(p) = 


M-fg /2 + p 0 — m — p 2 /2m + *0 


( 6 ) 


n (k) 


as 

47T 


C A n g (k) +T R niU q (k) , 


( 10 ) 


where Ca is the Casimir operator of the adjoint repre¬ 
sentation and T r is the trace normalization of the funda¬ 
mental representaion. We obtain the gluon contribution: 


n„,q = -il + r(£)e ™.(|l)‘ 

2 T(1 — e)T(| — e) 2 (5 — 4 e) F(2 — e ) 2 
r(| — 2e) + (l-p^r( 4 - 2 e) 

+ e {1 ~ P 2 (3 — 2 e)}/C(e, p) 


( 11 ) 


V [3 log(l —p 2 ) - 12 log 2 + 6 - 3 f{ P )] + 0(e), 

( 12 ) 
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(a) (b) (c) (d) 


/1C\ 

expansion of SE^ ot ) in as, it suffices to replace I(e, A) 
by its asymptotic expansion in A given by 

/(e, A) = /| 0 ) (e) + A • l£\e) + A 1_2e • J<& (e) + • • • . 

(18) 


FIG. 2: Feynman diagrams for the 1-loop Coulomb vacuum 
polarization: (a) massless quark loop, (b) transverse-Coulomb 
gluon (gr-gc) loop, and (c)(d) transverse gluon loops. 


where 7 F = 0.5772... is the Euler constant, p = fc 0 /|fc|, 
and 


K.(e,p) 

f(p) = 


= dx 
Jo J 

w 2 r 


dy 


1—e - —e 

x y 2 


[(1 - xy) - (1 - x) p 2 } 1+ ' 


; ( 13 ) 


Li 2 - 


P 

— 2 log 2 . 


1 ~P 
1 +P 


■log' 


( 1 + P 
\1-P 


7 r 2_l 
12 

(14) 


The analytic expression @ is a new result. Q On the 
other hand, the contribution from a massless quark is 
gauge independent and well known: 


n q (k) 






V 8 T (2 — e ) 2 
) T(4 — 2e) 

+ 0(e). 


(15) 

(16) 


For fcg > |fc| 2 the vacuum polarization II(fc) acquires 
an imaginary part; in this case, the usual +?'0 prescrip¬ 
tion, k,Q —> fcg + iO, is understood. If we take the in¬ 
stantaneous limit fcg —> 0 of the Coulomb propagator 
— C F (47ras/|fc| 2 ) [1 + II(fc)], the 1-loop static QCD po¬ 
tential in momentum space Jij is correctly reproduced. It 
is ensured by the Ward identity in Coulomb gauge [0 . 

Next we evaluate the integral One may easily ver¬ 
ify that the integral is finite both in ultraviolet and in¬ 
frared regions. Using the expressions (0): (11) and Jl5| ) 
for n(fc), we may take advantage of the functions J(e) 
and /(e, A) analyzed in [ju|. These functions are defined 
similarly to SE^Jp, where II(fc) is replaced by |fc| _2e or 
by (-£ 2fc o + l^l 2 - *°) _e in Eq. (§): 

(\k\~ 2e ) is = - CP °f m {C F a s m)~ 2e J(e), 

((-Z 2 k 2 0 + \k\ 2 -iO)- £ ) ls 

= - ^ Cf0 f m ( C F asm)~ 2e /(e, £A). (17) 

Here, A = C F as/ 2 represents the ratio of the Coulomb 
binding energy and the Bohr scale. Then SEp 0t can be 
expressed in terms of J(e, A) and J(e) plus an integral 

over x and y of terms involving J(e, y^f^A) and its 
derivative with respect to A. In order to obtain the series 


Only those terms relevant to our calculation are written 
explicitly. I Y and 7^2 were given in | pj|| , while I Y is 
needed additionally: 

4 1 ) (e)=4^- 3 / 2 r(l- e )r(i + e). (19) 


The integral over x and y can be evaluated analytically. 
We obtain 


V 7T 


(^) 2 (c A C F 4 2) + T Rni C F 4 2 )) 


+ 0(a%) 


20 ) 


with 



p A 

97 

+ 18’ 

( 21 ) 

C F asm J 

-l iog i 

( p N 

\ 22 

( 22 ) 

/ C F asm y 

/ 9 ’ 


d {2) = ^y - . 2 )log(C^ s ) - 10 + y + 7^3,(23) 

(24) 


4 2) = -g log(C F as) + 4. 


The 0(a 3 s m) terms [ 1 d Y and d^] stem from the static 
QCD potential and are well known. The 0(a A s m) terms 
[gY and 4Y are the new results. We made a cross 
check of our results through an independent computation 
using a different integral representation of n(7) and a 
slightly generalized version of the function /(e, A). We 
also checked our results by calculating the momentum- 
space integral (jjij) numerically. 

One may use the threshold expansion technique [0] in 
order to clarify from which kinematical regions individual 
perturbative corrections originate. We confirmed that, 
with respect to the quark-loop contribution, the potential 
region (ho ~ a 2 s m, |fc| ~ asm) accounts for ; the soft 
region induces —(8/3) log[/z/ {2C F asm)\ plus a constant; 
the ultra-soft region induces +(8/3) log[/Lt/(2C|.a|m)] 

plus a constant; the latter two contributions add up to 
( 2 ) 

produce dq '. Similar separate contributions can be iden¬ 
tified also for the gluonic contribution. 

( 2 ) 

As stated, the gluonic contribution d g y was computed 
numerically in ||. Comparing our result with the nu¬ 
merical results listed in TABLE I of that paper, we find 
that their results for the IS state are smaller than ours 
by about factor 6 , consistently for all the listed values of 
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C A C F d 

3.5 - 2.1 log(CVas) 

TrtiiCf d^ 

10.7-7.1 log(CVas) 

AQCDpot 

153.6 + 119.1 £ + 52.1 i 2 

A 1/r 2 

39.5 

Agreit (J = 1) 

-0.4 

-A]3reit (*/ = 0) 

23.0 


TABLE I: A numerical comparison of our results with the 
previously known corrections at 0{a%m). We take ni = 4. 
i = log \n/{CFasm)\. 

as- (Note that as in jo] should read C F as in contem¬ 
porary notations.) We have checked that Eq. (5.13) of 
that paper is correct. Therefore, we suspect that some 
error has occurred in transforming this equation to the 
final results given in the table of that paper. We are un¬ 
able to locate the error more precisely, since no details 
are provided for this part of the computation. 

Let us compare our results with the other 0{a‘gm) 
corrections: 

SE 1S = - j {C F as) 2 m x 

X (^d.QCDpot “b A\ / r 2 T ^Breit)- (25) 

We separate the corrections into 3 parts: AQCDpot de¬ 
notes the correction originating from the static QCD po¬ 
tential; Ai / r 2 denotes the correction originating from the 
—CAC F a 2 s /(2mr 2 ) potential; Aereit denotes the correc¬ 
tion originating from the Breit Hamiltonian. A numerical 


comparison is given in Table | for ni = 4 and for the stan¬ 
dard values of the color factors C F = 4/3, Ca = 3 and 
Tr = 1/2. We see that the corrections CUCf c?/ and 
T F niC F dq ; are not negligible as compared to A 1 / r 2 or 
A-Breit, whereas AQCDpot is an order of magnitude larger 
than the other corrections (for a typical choice of /r) due 
to an enhancement by O(Aqcd) renormalon p3| . 

The reason why the corrections calculated here have 
been overlooked in the effective theories appears to be 
two-folcl: (1) The present power counting schemes can¬ 
not specify all diagrams (and kinematical regions) which 
contribute to a given order of 1/c expansion. (2) The 
effective theories have been matched to on-shell QQ am¬ 
plitudes of perturbative QCD; generally the matching 
should be performed with off-shell QQ amplitudes or in¬ 
cluding on-shell amplitudes with additional gluons in ex¬ 
ternal lines [Q . We consider that we have not yet under¬ 
stood well higher-order corrections in the non-relativistic 
boundstate theory, especially where the soft and ultra- 
soft contributions are involved. 

Our results given here also apply to the QED bound- 
states composed of heavy particles such as the 
boundstate, after trivial replacements of the color fac¬ 
tors and the coupling constant. 

Acknowledgments 

Y. K. was supported by the Japan Society for the Pro¬ 
motion of Science. 


[1] A. Hoang, M. Smith, T. Stelzer and S. Willenbrock, Phys. 
Rev. D59, 114014 (1999); M. Beneke, Phys. Lett. B434, 
115 (1998). 

[2] K. Melnikov and A. Yelkhovsky, Phys. Rev. D59, 114009 
(1999); A. Hoang, Phys. Rev. D61, 034005 (200 0). 


[3] S. Komine and M. Yamaguchi, hep-ph/0110032 


381 (2001); 

heD-Dh/0108084 

f5l Y. Sumino. 

hep-ph/0104259 

ino, 

hep-ph/0109122 


Phys. Lett. B513, 
S. Recksiegel and Y. Sum- 


[6] A. Duncan, Phys. Rev. D13, 2867 (1976). 

[7] W. Fischler, Nucl. Phys. B129, 157 (1977). A. Billoire, 
Phys. Lett. B92, 343 (1980). 

[8] M. Peter, Phys. Rev. Lett. 78, 602 (1997); Nucl.Phys. 
B501, 471 (1997); Y. Schroder, Phys. Lett. B447, 321 
(1999); M. Melles, Phys. Rev. D62, 074019 (2000). 

[9] S. N. Gupta and S. F. Radford, Phys. Rev. D24, 2309 
(1981); Phys. Rev. D25, 3430 (1982). 

[10] S. Titard and F. J. Yndurain, Phys. Rev. D49, 6007 
(1994); Phys. Rev. D51, 6348 (1995); A. Pineda and F. 
J. Yndurain, Phys. Rev. D58, 094022 (1998); Phys. Rev. 
D61, 077501 (2000). 

[11] W. E. Caswell and G. P. Lepage, Phys. Lett. B167, 437 
(1986); B. A. Thacker and G. P. Lepage, Phys. Rev.D43, 
196 (1991); G. P. Lepage, L. Magnea, U. Magnea, C. 
Nakleh and K. Hornbostel, Phys. Rev. D46, 4052 (1992). 

[12] A. Pineda and J. Soto, Nucl. Phys. Proc. Suppl. 64, 428 


(1998). 

[13] M. Beneke and V. Smirnov, Nucl. Phys. B522, 321 
(1998). 

[14] N. Brambilla, A. Pineda, J. Soto and A. Vairo, Phys. 
Lett. B470, 215 (1999); B. Kniehl and A. Penin, Nucl. 
Phys. B563, 200 (1999); Nucl. Phys. B577, 197 (2000). 

[15] Y. Kiyo and Y. Sumino, Phys. Lett. B496, 83 (2000). 
16 A. Hoang, A. Manohar and I. Stewart, Phys. Rev. D64, 

014033 (2001). 

[17] W. Kummer and W. Modritsch, Z. Phys. C66, 225 
(1995). 

[18] Y. Kiyo and Y. Sumino, hep-ph/0110277. 

[19] See e.g.: G. T. Bodwin and D. R. Yennie, Phys. Rept. 
43, 267 (1978). 

[20] S. Love, Ann. of Phys. 113, 153 (1978). 

[21] D. Heckathorn, Nucl. Phys. B156, 328 (1979); G. Ad¬ 
kins, Phys. Rev. D27, 1814 (1983). 

[22] F. Feinberg, Phys. Rev. D17, 2659 (1978). 

[23] U. Aglietti and Z. Ligeti, Phys. Lett. B364, 75 (1995). 

[24] K. Hasebe and Y. Sumino, Phys. Rev. D61, 105001 

( 2000 ). 

[25] As for the main results A v M n i and A se A/„; of ||, the 
sign of the former correction is wrong; we discuss the 
error of the latter correction below. 

[26] A one-parameter integral representation is given in Q. 
















